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I ' We use quasi-simultaneous near-infrared (near-IR) and optical spectroscopy from 

2 I four observing runs to study the continuum around 1 /im in 23 well-known broad- 

emission line active galactic nuclei (AGN). We show that, after correcting the optical 
^ ' spectra for host galaxy light, the AGN continuum around this wavelength can be ap- 

proximated by the sum of mainly two emission components, a hot dust blackbody 
and an accretion disc. The accretion disc spectrum appears to dominate the flux at 
~ 1 ^m, which allows us to derive a relation for estimating AGN black hole masses 
^ , based on the near-IR virial product. This result also means that a near-IR reverber- 

ation programme can determine the AGN state independent of simultaneous optical 
spectroscopy. On average we derive hot dust blackbody temperatures of ~ 1400 K, a 
value close to the sublimation temperature of silicate dust grains, and relatively low 
hot dust covering factors of ~ 7%. Our preliminary variability studies indicate that in 
most sources the hot dust emission responds to changes in the accretion disc flux with 
the expected time lag, however, a few sources show a behaviour that can be attributed 
to dust destruction. 
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1 INTRODUCTION rever beration mapping studies (see, e.g., review bv lPetersonI 

, , .... . „ . , . Il993l l. We have started to exten d these stud i es to near- 

ine broad emission line region (BbLH) oi active galactic . „ , , T.r>\ ^ t IT 71 — I — Tl JnorvcJ i 

, . . r 1 )■ r , ■ infrared (near-lH) wavelengths. In lLandt et al.l (I2008I . here- 

nuclei I AGN I IS one oi the most direct tracers ot the imme- r, t-. t\ .i- i .i ^ r 

^ ^ , , , , , after Paper 1) we outlined the rationale ot our programme, 

diate environment of supermassive black holes. However, de- ^1^11 r J2 ^ ^ 11 

presented the observations of the farst three epochs and ad- 

spite decades of intensive optical and ultraviolet (UV I spec- , , , . n r,i • , . ■ i- ,1 

^ ' dressed brietfy some of the important issues regarding the 

trophotometric studies its geometry and kinematics remain , . , . , , , .... 

, , I-— — -11 -1 physics ot the most prominent broad emission lines. Here we 

lU-defaned (see, e.g.. review bv lbulentic et al.ll200q i. Our cur- ^ r ii 1 r i- j ■ ^ 



present the fourth epoch of observation and investigate the 
continuum around the 1 nm inflection point. 



rent, limited knowledge of its physical condition and scale 

was gained primarily through th e application o f photoionisa- 

tion models (see, e.g., review by |Ferland|[20o3 ) and through The AGN spectral continuum region around the rest- 

frame wavelength of ~ 1 /im is believed to sample simultane- 
ously two impor tant emission componen t s, namely, th e ac- 
E-mail: hlandt@unimelb.edu.au cretion disc fe.g.. lMalkan fc Sarge"iitlll982l : [Malkanlll983l') and 



t Visiting Astronomer at the Infrared Telescope Facility, which is the hottest part of the putativ e dusty torus fe.g.. iBarvainT^ 

operated by the University of Hawaii under Cooperative Agree- |^. [Neugebauer etlol fl987l ) . However, although it is as- 



ment no. NNX-08AE38A with the National Aeronautics and 



^ . , . . . „ . . T^. T,, A sumed to be understood, it has not yet been sampled spec- 

bpace Administration, bcience Mission Directorate. Planetary As- t-, i 1 1 

tronomy Program. ' troscopically m its entirety. By probi ng the long- wavelength 

t Current address: Department of Physics and Astronomy, Geor- end o f the accretion disc spectrum ([Kishimoto et al.||2005|, 

gia State University, 709 One Park Place South, Atlanta, GA l2008l ) such an investigation has the potential to solve the 

30303, USA discrepancy often found between theoretical models and ob- 
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servations (see, e.g., review bv lKoratkar fc Blaes|[l999l ) . Fur- 
thermore, since the accretion disc most hkely iUuminates 
directly the inner parts of the dust structure, monitoring 
the change in spectral flux and slope of these two emission 
components relative to each other can constrain the loca- 
tion and geometry of the obscurer. So far, only few sources 
have been observed in d ust reverberation progra mmes (e.g. . 
Glass" 1992'; 'Nelson" 1996; O knvanskii e t al. 1999: Glass" 2004 
Minezaki et al.. .2004; .Suganuma et al.. ,2006, ; .Koshida et al.l 
20091 ). 

The paper is organised as follows. In Section 2 we briefly 
introduce the sample and discuss the observations. In Sec- 
tion 3 we derive pure AGN continuum spectral energy distri- 
butions (SEDs) , based on which we constrain the individual 
continuum components (Section 4). The variability of these 
components is disussed in Section 5. Finally, in Section 6 
we summarize our main results and present our conclusions. 
Throughout this paper we have assumed cosmological pa- 
rameters Ho — 70 km s~^ Mpc~^, Q,m = 0.3, and Qa ~ 0.7. 



2 THE OBSERVATIONS 

The target selection, observational strategy, and data re- 
duction procedures have been described in detail in Paper 
I. In short, we obtained for a sample of 23 well-known rel- 
atively nearby [z < 0.3) and bright (J < 14 mag) broad- 
emission line AGN during four observing runs contempora- 
neous (within two months) near-IR and optical spectroscopy. 
The observations were carried out between 2004 May and 
2007 January with a single object being typically observed 
twice within this period. 

In the near-IR we used the SpeX spectrograph 
l|Ravner et all 120031 ') at the NASA Infrared Telescope Fa- 
cility (IRTF), a 3 m telescope on Mauna Kea, Hawai'i. We 
chose the short cross-dispersed mode (SXD, 0.8 — 2.4 ^m) 
and a short slit of 0.8 x 15". The op tical spectra were ob - 
tained with the FAST spectrograph l|Fabricant et aLlFlGOSl ) 
at the Tillinghast 1.5 m telescope on Mt. Hopkins, Arizona, 
using the 300 1/mm grating and a 3" long-slit, resulting in 
a wavelength coverage of ~ 3720 - 7515 A. 

Paper I presented details for the first three epochs (2004 
May, 2006 January, and 2006 June). Similarly, we give in Ta- 
bles [1] and [2] the details for the fourth IRTF epoch carried 
out on 2007 January 24-26. These nights were mostly pho- 
tometric with a seeing in the range of ~ 0.7 — 1". 



3 AGN SPECTRAL ENERGY 
DISTRIBUTIONS 

Our aim is to study the spectral shape of the AGN contin- 
uum over the entire observed frequency range and how it 
changes as the AGN flux varies. For this purpose we need 
to isolate the spectrum of the pure AGN and ensure an ac- 
ceptable relative photometry between epochs. In practice 
this means that we need to subtract the flux of the under- 
lying host galaxy and correct the individual spectra for flux 
calibration errors (e.g., varying weather conditions, aperture 
and seeing effects). In the absence of short-term variability, 
these steps should yield an alignment between the AGN op- 
tical and near-IR spectral parts. 
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Figure 1. Spectral aperture versus enclosed host galaxy luminos- 
ity at rest-frame 5100 A as derived from HST images. Based on 
the observed correlation (solid line) we have estimated the host 
galaxy flux in the spectra of sources without useful HST images. 
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Figure 2. Histogram of the ratio between the optical correction 
factors applied to the host galaxy and the total flux. 



In principle, accurate (within a few per cent) relative 
photometry of the spectra can be achieved a posteriori by us- 
ing the fact that the narrow emission line fluxes o f AGN are 
constant over several years l|Petersonlll98i , 1 19931 ) . However, 
this requires that the observations were performed under 
similar conditions, i.e., using the same aperture and slit po- 
sition angle and having the same seeing. In general, these re- 
quirements were not met for our optical observations, which 
were obtained in service-mode under widely varying weather 
conditions. On the other hand, the near-IR spectra were ob- 
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Table 1. IRTF Journal of Observations for 4th Epoch 



Object Name observation exposure airmass continuum S/N telluric standard star 





date 






[sec] 






J 


H 


K 


name 


distance 
[deg] 


airmass 


(1) 


(2) 






(3) 




(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Mrk 335 


2007 


Jan 


25 


32x120 


1 


.431 


33 


110 


272 


HD1160 


16.1 


1.977 


Mrk 590 


2007 


Jan 


24 


58x120 


1 


.214 


37 


122 


143 


HD13936 


10.2 


1.697 


Ark 120 


2007 


Jan 


26 


48x120 


1 


.162 


187 


226 


436 


HD34317 


1.8 


1.054 


Mrk 79 


2007 


Jan 


25 


48x120 


1 


.334 


97 


202 


448 


HD45105 


12.6 


1.128 


PC 0844+349 


2007 


Jan 


24 


48x120 


1 


.138 


99 


170 


259 


HD71906 


4.2 


1.271 


Mrk 110 


2007 


Jan 


26 


48x120 


1 


.449 


61 


148 


452 


HD71906 


17.8 


1.368 


NGC 3227 


2007 


Jan 


25 


16x120 


1 


.339 


57 


94 


146 


HD89239 


7.6 


1.122 


NGC 4151 


2007 


Jan 


24 


8x120 


1 


.446 


15 


54 


113 


HD109615 


4.9 


1.398 


3C 273 


2007 


Jan 


25 


40x120 


1 


.060 


123 


343 


347 


HD109309 


11.6 


1.199 


HE 1228+013 


2007 


Jan 


25 


48x120 


1 


.370 


40 


101 


274 


HD109309 


10.7 


1.242 


NGC 4593 


2007 


Jan 


24 


8x120 


1 


.446 


43 


123 


191 


HD112304 


10.8 


1.540 


NGC 5548 


2007 


Jan 


24 


16x120 


1 


.319 


57 


73 


145 


HD131951 


14.0 


1.299 


Mrk 817 


2007 


Jan 


26 


40x120 


1 


.509 


98 


224 


323 


HD121409 


7.8 


1.480 



The columns are: (1) object name; (2) observation date; (3) exposure time; (4) average airmass; S/N in the 
continuum over ~ 100 A measured at the central wavelength of the (5) J, (6) H, and (7) K band; for the star 
used to correct for telluric absorption (8) name, (9) distance from the source, and (10) average airmass. 



Table 2. Tillinghast Journal of Observations for 4th Epoch 



Object Name 


IRTF 




observation 


exposure 


airmass 


cloud 




observation 


date 






[sec] 






condition 




date 




















(1) 


(2) 






(3) 






(4) 


(5) 


(6) 


Mrk 335 


2007 


Jan 


25 


2007 


Jan 


24 


2x300 


1, 


,38 


clear 


Mrk 590 


2007 


Jan 


24 


2007 


Jan 


24 


2x720 


1, 


,21 


clear 


Ark 120 


2007 


Jan 


26 


2007 


Feb 


08 


2x480 


1, 


,18 


cirrus 


Mrk 79 


2007 


Jan 


25 


2007 


Jan 


22 


2x480 


1, 


,12 


clear 


PG 0844+349 


2007 


Jan 


24 


2007 


Jan 


22 


2x480 


1, 


,03 


clear 


Mrk 110 


2007 


Jan 


26 


2007 


Jan 


22 


2x420 


1, 


,08 


clear 


NGC 3227 


2007 


Jan 


25 


2007 


Jan 


22 


2x180 


1, 


,02 


clear 


NGC 4151 


2007 


Jan 


24 


2007 


Feb 


12 


2x 90 


1, 


,07 


cloudy 


3C273 


2007 


Jan 


25 


2007 


Feb 


12 


2x300 


1, 


,25 


cloudy 


HE 1228+013 


2007 


Jan 


25 


2007 


Feb 


12 


2x480 


1, 


,30 


cloudy 


NGC 4593 


2007 


Jan 


24 


2007 


Feb 


17 


2x270 


1, 


,27 


cirrus 


NGC 5548 


2007 


Jan 


24 


2007 


Feb 


17 


2x270 


1, 


,02 


cirrus 


Mrk 817 


2007 


Jan 


26 


2007 


Fob 


17 


2x360 


1, 


,13 


cirrus 



The columns are: (1) object name; (2) IRTF observation date (reproduced from Ta- 
ble[T]l; (3) Tillinghast UT observation date; (4) exposure time; (5) average airmass; 
and (6) cloud condition. 



tained uniformly during four separate epochs, with one of 
them (2006 June) having photometric weather conditions 
throughout. In addition, the much better seeing achieveable 
in the near-IR relative to the optical band allowed us to use 
a relatively small slit for these spectra, which minimized the 
contamination by host galaxy light. 

Given these considerations, we chose the following ap- 
proach to obtain pure AGN continuum SEDs with accept- 
able relative photometry. Firstly, we performed relative pho- 
tometry of the near-IR spectra using the flux of the strong 
narrow emission line [S III] A9531 with the epoch observed 
under the best weather conditions (usually 2006 June or 
2007 January) as a reference. We have measured the emis- 
sion line fluxes as described in Paper I and list them in Table 
[3] We note that [S ill] A9531 is not observed for the source 
HE 1228+013, and therefore, no adjustment of the near- 



IR spectra was made in this case. Secondly, we subtracted 
the host galaxy flux from the (adjusted) near-IR spectra. 
Finally, we corrected the flux scale of the optical spectra 
by requiring that after the subtraction of the (optical) host 
galaxy contribution the optical and near-IR AGN spectral 
parts align. In Fig. lAll we show the original data (left pan- 
els) and the data after host galaxy subtraction and spectral 
alignment were applied (right panels). 

We have estimated the host galaxy contribution in the 
apertures of both the near-IR and optical s pectra using 
the Hubble Space Telesco p e (HS T) images of iBentz et al.l 
(2006E) and iBentz et al] (|2009l ) and following their ap- 
proach. The observed HST fluxes were transformed to a 
rest-frame wavelength of 5100 A by applying a color cor- 
rection factor based on the model bulge galaxy template of 
iKinnev et all l|l996l ) and were corrected for Galactic extinc- 
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Figure 3. Histograms of the ratio between tlie line fluxes corre- 
sponding to the best-weather IRTF run and another epoch. The 
top, middle and bottom panels show the narrow emission lines 
[S III] A9531, [O III] A5007, and [O III] A5007 with the optical 
correction factor listed in Table |4] applied, respectively. 




2 3 
opt. corr. factor 

Figure 4. The ratio between the [O III] A5007 line fluxes from 
optical reverberation programmes and our observations versus the 
correction factor applied to the optical spectrum. The locus of 
equality is shown as the solid line. 



tion usin g A\ values derived f rom th e hydrogen column den- 
sities of Dickey fc LockmanI l|l990t ). The unabsorbed rest- 
frame 5100 A fluxes were then use d to scale the galaxy tem- 
plate from iMannucci et al.l (|200ll ') of the appropriate Hub- 
ble type, which was subsequently subtracted from the (rest- 
frame) spectrum. Details of our host galaxy flux estimates 
are given in Tables |3] and |4] for the near-IR and optical spec- 



tra, respectively. As expected, the relatively small aperture 
used in the near-IR includes significantly less host galaxy 
flux than the optical aperture, on average a factor of ~ 3. 



For 8/23 objects we did not have suitable HST images. 
In these cases we have estimated the host galaxy contri- 
bution based on the linear correlation present between the 
logarithms of the aperture and the enclosed host galaxy lu- 
minosity for the sources with HST images (Fig.[T|. The ob- 
served correlation is logLhost = 0. 73(±0. 04) dog (aperture) -f- 
39.05(±0.04), where Lhost is the host galaxy luminosity at 
rest-frame 5100 A (in erg s~^ A~^) and aperture is the spec- 
tral aperture (in kpc^). Additionally, we could not find in- 
formation on the host galaxy type for four sources. In these 
cases we have assumed elliptical hosts for the two high- 
redshift sources IRAS 1750-1-508 and PDS 456 and spiral 
(SO) hosts for the two low-redshift sources H 2106—099 and 
H 1934-063. 



The frequency gap between our near-IR and optical 
spectra is large enough (log v ~ 0.05) that we need to as- 
sume an overall spectral shape in order to judge the two 
parts as being 'aligned'. Therefore, since the AGN contin- 
uum blueward of ~ 1 fj,m is generally thought to be emitted 
by the accretion disc, we have considered this component as 
the model and calculated its spectrum as detailed in Sec- 
tion |4]T] Assuming an accretion disc spectrum, excess host 
galaxy contribution to the optical spectrum (but not to the 
near-IR spectrum) was apparent in the majority of our low 
redshift [z <. 0.1) sources, requiring also a host galaxy flux 
correction factor. However, as Fig. [5] shows, in most cases 
the host correction factors are similar to the correction fac- 
tors applied to the total optical spectra, and in only a few 
sources they indicate a real underestimation (by factors of 
<3). 



The relative photometry in the near-IR required flux 
correction factors of < 2, with the exception of Mrk 590 (fac- 
tor of ~ 3), and most spectra had flux differences relative to 
the best- weather epoch of no more than ^ 20% — 30% (Fig. 
[31 top panel). The correction factors for the optical spectra 
spanned a larger range of ~ 0.6 — 3.8. The narrow emission 
line [O III] A5007 is usually used for relative photometry in 
the optical and we list in Table U also its flux. As Fig. [S] 
shows, whereas our approach did not succeed in bringing 
also the [O III] A5007 line fluxes of the different epochs into 
alignment (bottom panel), it did reduce the scatter of the 
original flux ratio distribution (middle panel). Finally, for a 
subsample of 12 sources we compare in Fig.Uthe optical cor- 
rection factors with the ratio between the [O ill] A5007 line 



fluxes observed by reverberation programmes i Korista et al 



19951 : iPeterson et al. l ll998l : lBentz et al. I l2006al : 



Dennev et al 



20101) and us. Fig.|4]indicates that the optical correction fac- 



tors are mostly overestimates relative to photometric condi- 
tions. This result could be partly due to our near-IR flux- 
calibration, which was based on the telluric standard stars. 
Since the target integration times were considerably larger 
than those of the stars, the seeing would have varied sub- 
stantially for the former, leading overall to an overestimate 
of the extraction aperture and so of the flux scale. 
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Table 3. Estimates of the Near-IR Host Galaxy Contribution 



Object Name 


z 


A{l+z)5100 


host 


ref. 


IRTF 


[S III] A9531 


near-IR spectrum 


HST flux 


HST flux 






[mag] 


type 




run 


Flux 


aperture 


PA 


(5580 A) 
[erg/s/cm^/A] 


((1 + 2)5100 A) 














[erg/s/cm^] 


[arcsec] 


[deg] 


[erg/s/cm^/A] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




(8) 


(9) 


(10) 


(11) 




u.ouu 




? 




9004 A/Tn-ir 


Q 88F 1 ^ 


0.8 


X 


3.0 







4 Q'?P 1 7 
















0.8 


X 


4.0 







R ORP 1 7 


J.J_ i O --^ i 1 U O 


0.297 


0.061 




F04 


0004 IVFav 


1.25E— 14 


0.8 


X 


3.2 







5.29E— 17 




0.184 


0.941 


7 




9nnfi Tnn 


9 SI p,_i 

Zr. O ±i_J J-U 


0.8 


X 


4.0 







1 IQP—IR 

-L . -L yj_j — ±u 


3C 273 


0.158 


000 


E 


B09 


90nfi Tan 


1.13E— 14 


0.8 


X 


5.4 


207 


9 Rl F_l (\ 


9 c;'^P_l R 

Z . UtJJ-J — -LU 












9nnfi Tnn 




0.8 


X 


4.0 







9 99P_1 R 












9007 Tan 


1 48P,— 14 


0.8 


X 


5.4 





9 c;OP— 1 R 


9 44P— 1 R 

Zr.'rt'iJJj — -LU 


A/Trlr 87fi 


0.129 


0.005 


E 


B09 


9nnfi Tnn 


1.17E— 14 


0.8 


X 


4.4 





A (SOP— 1 R 


4 '^4P— 1 R 

1. Urli-J — -LU 


HP, 1 9984-01 


117 


000 


E 


L07 


9nnfi Tnn 

^\J\J\J iJ Llll 

9007 Tan 




0.8 
0.8 


X 
X 


3.4 
5.4 








1.75E— 16 

9 41 P— 1 R 


J. VJ VJ O rt'± ij 'rt 


0.064 


0.031 


Sa 


B09 


2006 Tan 

i*VJL/VJ U till 


4.57E— 15 


0.8 


X 


5.0 


270 


c; 27F— 1 6 


nnF— 1 6 

U.UWl-J -LU 












9nnT Tan 

^LJLF 1 (J cLll 


7.10E— 15 


0.8 


X 


4.6 





c; 'ifiP—i R 

tj. OUH/ — -LU 


c; 08P— 1 R 

U. UOH/ — -LU 


A/Trlr 1 1 D 


035 


000 


Sa 


P07 


900fi Tan 


1 70F,— 14 


0.8 


X 


5.4 


254 


9 4QP— 1 fi 


9 1 SiP_1 R 

Zi . J-OJ-J — -LU 












9007 Tan 
.Ld\}yy 1 1 1 cLii 


1 q7F,_14 


0.8 


X 


5.0 





9 47P— 1 fi 


9 1 RP_1 R 

Zi . J-UJ-J — -LU 


A/Trlr '^OQ 
iV±l rv oyjij 


0.034 


0.083 


SO /a 


P07 


9004 A/lav 




0.8 


X 


2.8 





O.OOi-J — i-U 


SI P_1 R 












9nnfi Tnn 

^LJLFU (J Llll 


■4:. .^!jl_J A.'-t 


0.8 


X 


4.0 







4 77P_1 R 


Ark 1 2(1 


0.033 


0.570 


Sb 


P07 


2006 Tan 


2.61E— 14 


0.8 


X 


4.0 


327 


4.10E— 15 


tj. O 1 1-J -LU 












9007 Tan 

^\J\J I 1 1 cLll 


9 98P,— 1 4 


0.8 


X 


5.0 





4.11E— 15 


O.OOi-J -LU 


Mrk 81 7 


0.031 


0.000 


Sa 


P07 


9004 Mav 


2 '^SF,— 14 


0.8 


X 


3.6 





2 26F— 1 6 


1 QRP— 1R 

-L.C/Ul-J -LU 












900fi Tnn 

^ULFU (J Llll 


9 98P,— 1 4 


0.8 


X 


4.0 





9 c;i P_1 R 

Zr. u ± jjj — J_U 


9 1 SP_1 R 

Zr. J- OH/ — -LU 












2007 Tan 


2.70E— 14 


0.8 


X 


5.0 





2 Q^F,— 1 6 

Zj.C/t^l— J J-U 


2 ^7F— 1 6 

Zj • U f l-J -L U 


A/TrW 9Qn 


030 


000 


SO 


P07 


9004 A/Tnv 


9 08F,— 1 4 


0.8 


X 


3.0 







4 71 P_1 R 

4:. 1 ±J_LJ -LU 












900fi Tnn 

^\J\J\J iJ Llll 


^'iF,_1 4 

O. OOi-J — J-rt 


0.8 


X 


4.0 







^ SOP— 1 R 

U. OUH/ — XU 




0.027 


0.161 


7 




9004 Mav 


2 '^^F— 14 


0.8 


X 


3.6 







5.74E— 16 












900fi Tnn 

^\JWLJ iJ Llll 


2 25F— 1 4 


0.8 


X 


4.0 







R 9QP— 1 fi 


Mrk S^"! 


0.026 


0.077 


so/a 


P07 


2006 Tan 


1.31E— 14 


0.8 


X 


4.4 


86 


>S Q7F— 1 6 

<J 'iJ i 1-J -LU 


5.44E— 16 












2007 Tan 

^\ } yj 1 «j dii 


1.28E— 14 


0.8 


X 


4.2 





fi 21 F— 1 6 

U*Z*-l1._J -lU 


U. U 1 l-J -LU 


A/Trlr '^QD 


0.026 


0.000 


SO 


P07 


900fi Tan 


1 10P,— 14 


0.8 


X 


5.0 


353 


1.31E— 15 


l.llE— 15 












9007 Tan 
^xjyj 1 (J dii 




0.8 
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Figure 5. Left panels: Rest-frame AGN (host galaxy-subtracted) spectral energy distributions (from right panels of Fig. lAljl normalized 
at 1 lira (vertical dashed line) to the luminosity of the lowest-flux epoch. The four IRTF observing runs are shown in colour: 2004 May 
(black), 2006 January (green), 2006 June (red), and 2007 January (blue). The dotted line indicates the accretion disc spectrum that 
approximates best the continuum blueward of ~ 1 /^m. The horizontal dashed line marks the minimum in the integrated luminosity 
of the lowest-flux epoch. Wavelength units in fira are labeled on the top axis. Right panels: As in left panels with the accretion disc 
spectrum subtracted. The dotted line shows the blackbody spectrum fitted to the continuum redward of ~ 1 ^ira. 
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Table 4. Estimates of the Optical Host Galaxy Contribution 



Object Name z A(i_,.2)gioo IRTF [O III] A5007 opt. optical spectrum HST flux HST flux host 
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The columns are: (1) object name; (2) redshift from the NASA/IPAC Extragalactic Database (NED); (3) Galactic extinction at rest-frame 5100 A; 
(4) IRTF run; (5) observed flux of the narrow-line [O III] A5007; (6) flux scale correction factor for optical spectrum; (7) optical extraction aperture; 
(8) optical slit position angle, where PA=90° correspo nds to E-W orientat i on; (9 ) observed host galaxy flux at 5580 A (except for Mrk 509 at 
5483 A) in optical aperture derived from HST images of lBentz et al. I l l2006bl . [200g) ; (10) observed host galaxy flux at rest-frame 5100 A in optical 
aperture, corrected for Galactic extinction using the values in column (3); and (11) host galaxy flux correction factor. 



4 THE AGN CONTINUUM COMPONENTS 

Our observations cover the continuum rest-frame frequency 
range oi u ^ 10^'* — lO'^^ Hz, which is believed to sam- 
ple two import ant emission components , namely, the a ccre- 
tion disc (e.g., iMalkan fc SargentI [l98^ : iMalkanI [l983l ) and 



the hottest part of the putativ e dusty torus (e.g.. iBarvainisI 
1 19871 : iNeugebauer et al.lll987l ). These two components are 
predicted to have opposite spectral behaviours, namely, the 
accretion disc and the hot dust emission rising and falling 
towards longer frequencies (in a logarithmic vf^ versus u 
plot), respectively, resulting in an inflection point at the lo- 
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Figure 5 - continued 



cation where they meet (at ^ 1 um; e .g.. IWard et al. 



1987 



200e 



Carleton et al. Ill987|; lElvis et al.lll994l : iGlikman et al, 
Riff cl ct al. 2006) 



Our present knowledge of this special continuum region 
is based mainly on optical spectroscopy combined with near- 
IR photometry, with the two frequency ranges usually not 
observed contemporaneously. Our observations now give us 
for the first time the opportunity to define the separate con- 



tinuum components at spectroscopic rather than photomet- 
ric precision. This permits us to investigate their properties 
in greater detail than was possible before. In particular, we 
are interested to understand if a near-IR reverberation map- 
ping campaign generally requires simultaneous optical spec- 
troscopy to determine the state of the ionising fiux or if the 
latter can also be derived from the continuum behaviour at 
longer wavelengths. 
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Figure 5 - continued 



For the purpose of the following analysis we have nor- 
malised all SEDs at ~ 1 /im (Fig. The SEDs of two 
sources, namely, Mrk 590 and NGC 3227, are strongly af- 
fected by host galaxy light in both their near-IR and optical 
spectral parts (Fig. [6| and are not considered further. 



4.1 The accretion disc 

The AGN continuum blueward of ~ 1 ^m is believed to be 
emitted by the accretion disc, which is considered the main 
source of ionising flux producing the broad emission lines. 
In order to test if this component is indeed seen in our data, 
we have calculated accretion disc spectra. We have assumed 
a steady geometrically thin, optically thick accretion disc, in 
which case the emitted flux is independent of viscosity and 
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Figure 5 - continued 



each element of the disc face radiates roughly as a black- 
body with a characteristic temperature depending only on 
the mass of the black hole, Afan, the accreti on rate, M, and 
the radius of the i nnermost stable orbit fe.g.. |Peterson|[l997l : 
iFrank et al.ll200^ . We have adopted the Schwarzschild ge- 
ometry (non-rotating black hole) and for this the innermost 
stable orbit is at rin = 6 • rg, where r-g is the gravitational 
radius defined as rg = GMbh/c^, with G the gravitational 



constant and c the speed of light. Furthermore, we have as- 
sumed that the disc is viewed face-on. 

The accretion disc spectrum is fully constrained by 
the two quantities, mass and accretion rate of the black 
hole. Two thirds of our sample have black hole masses de- 
rived from reverberation mapping campaigns. For the re- 
mainder (7 sources) we have estimated this quantity by 
applying the virial theorem Mbh v'^r/G, where v and 
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Figure 5 - continued 



r are the velocity and radial distance of an orbiting par- 
ticle, respectively, to the near-lR. Using the width of the 
Pa/3 broad component (denoted FWHMpa;3, published in 
Paper I) as a measure of v and the square-root of the con- 
tinuum luminosity at 1 /im (denoted J^Lipm) as a surro- 
gate for r, Fig. [7] shows that for the sources with rever- 
beration mapping results the black hole mass correlates 
with the near-IR virial product. The observed correlation 



is log A/bh = 0.84 ■ (2FWHMpa;3 + 0.5^Lip,„) - 16.58. The 
accretion rate can be obtained directly from an approxima- 
tion of the accretion disc spectrum to the data. We show 
our results in Fig. [S] left panels (dotted lines) and list the 
relevant values in Table [5] 

The calculated accretion disc spectrum approximates 
well the AGN continuum slope blueward of ~ 1 /xm in all 
high-redshift [z 0.1) sources, where the host galaxy con- 
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Figure 5 - continued 

tribution to the total flux is negligible at all wavelengths. disc emission can still dominate at the large wavelength of 
In particular the near-IR spectral part (~ 0.8 — 1 ^m, cov- ~ 1 ^m. This means that we observe the flux of the ionising 
ered by the IRTF spectrum) is well reproduced, which is component directly in the near-IR, which will allow us to 
free from major contaminating components such as, e.g., determine the AGN state without optical spectroscopy. 
strong Fe II emission and the 'small blue bump' (iGrandil 

ll982l : IWills et al.lll985l ) sometimes found in the optical spec- In the lower redshift sources, the host galaxy contri- 

tral part. However, most important for our future near-IR bution is expected to be minimal in the near-IR spectrum, 
reverberation programme is the result that the accretion since these observations were taken through a relatively nar- 
row slit, and in the bluer part of the optical spectrum, since 
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Figure 5 - continued 



at these frequencies the AGN emission lies well above that 
of the host galaxy. Indeed, in both these frequency regimes 
the calculated accretion disc spectrum is always a satisfac- 
tory approximation. However, excess host emission is often 
observed at the red end of the optical spectrum, which has 
the effect of flattening the spectrum (in a logarithmic uf,, 
versus v plot). Nevertheless, subtraction of sufficient host 



galaxy ffux leads also in this spectral range to an alignment 
of the data with the expected accretion disc emission. 



4.2 The hot dust 

The AGN continuum redward of 1 /xm is believed to be 
produced by the hottest part of a dusty torus, which sur- 
rounds the central ionising source and obscures the BELR 



14 H. Landt et al. 



Mrk 590 (z = 0.026) NGC 3227 (z = 0.004) 

2 1 0.8 0.6 0.4 2 1 0.8 0.6 0.4 




14 14.2 14.4 14.6 14.8 15 14 14.2 14.4 14.6 14.8 15 



log I'resi [Hz] 

Figure 6. Rest-frame spectral energy distributions for the two sources (Mrk 590 and NGC 3227) ttiat are strongly affected by tiost 
galaxy light. The spectra corresponding to the IRTF o bserving runs in 2006 J anuary (green) and 2007 January (blue) are normalized 
at 1 lira (dashed line). The host galaxy templates from lMannucci et al.l ll200ll ) are ovcrplotted (magenta). Wavelength units in /im are 
labeled on the top axis. 
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Figure 7. Black hole mass determined from reverberation cam- 
paigns versus the near-IR virial product between the width of the 
Pa/3 broad component and the continuum luminosity at 1 lira. 
Based on the observed correlation (solid line) we have estimated 
the black hole mass of sources without reverberation results. 



for lines of sight close to the plane of the accretion disc. As- 
suming this component is indeed seen in our data, we have 
subtracted from the total spectrum that of the accretion disc 
and have fitted to the result a blackbody spectrum. For this 
purpose we have used the C routine MPFIT (version 1.1; 
lMarkwardtll2009l ). which solves the least-squares problem 
with the Levenberg-Marquardt technique, and have fitted 
for the temperature and flux scaling. We have included in 
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Figure 8. Histogram of the temperature of the hot dust. 



the fit only the continuum part of the near-IR spectrum, 
i.e., we have excluded emission lines, and rebinned it to 
A log V — 0.01 Hz. We show our results in Fig.[5l right panels 
(dotted lines) and list the relevant values in Table [S] 

A hot blackbody spectrum appears to approximate 
well the near-IR AGN continuum in all our sources. 
In particular, the typical curvature of such a spectrum 
is evident in our spectra due to their relatively large 
wavelength coverage. This b e havio ur was noted also by 
iRodrfguez-Ardila fc Mazzalavl ()2006h . whose near-IR spec- 
trum of one broad-line AGN (Mrk 1239) covered an even 
larger wavelength range of 0.8 — 4.5 /im. In this respect. 
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Table 5. AGN Continuum Fit Parameters 
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The columns are: (1) obj ect na me; (2) black hole rn ass (in solar masses); (3) re ference for the black hole ma ss, where 
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we note that the blackbody curvature is not evident in 
the original spectra (without the accretion disc compo- 
nent subtracted), which resemble rather a single power-law. 
The resulting temperatures for the hot blackbody compo- 
nent are in the range of Thot ~ 1100 — 1700 K, which 
are typical values of the dust sublimation temperature for 
most astrophy sical grain compositions (« 1000 — 2000 K; 
ISalpeteilll977l ). The overall temperature distribution is rel- 
atively narrow and has a well-defined mean of (Thot) = 
1365 ±18 K (Fig. [S]). Similarly narrow temperature distribu- 



tions were obtained previously by AGN s tudies using low- 
and medium-resolution near-IR spectra l|Kobavashi et all 
ll993l : lRifrel et al.ll2009l ) and more recently by studies based 
on high-resol ution near-IR photometric and inte rferometric 
observations (|Kishimoto et al.ll2007l . |2009| . I2OI0I ') . 

This well-defined peak of the blackbody temperatures is 
close to the sublima tion temperature of silicate dust grains 
(|Kimura et al.|[2"002l ). More refractory dust types are known, 
notably carbonaceous dust, e.g. graphite, which can survive 
up to ~ 2000 K. The absence of such hot dust, except in one 
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Figure 9. Histogram of the hot dust covering factor C ~ 0.4 ■ 
(i/Lhot/'^i'acc), where fLhot and i^Lacc are the peak luminosities 
of the hot blackbody and accretion disc spectrum, respectively. 
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Figure 10. The hot dust covering factor versus the correction 
factor applied to the optical spectrum. 



object (NGC 5548), suggests that either carbonaceous dust 
is rare in AGN or that some other mechanism than dust sub- 
limation sets the maximum dust temperature. Dust forma- 
tion is critically dependent on the numerical ratio of carbon 
to oxygen atoms in th e parent gas; m ost pair up to form CO, 
and do not form dust (|Whittetli200^ 1. An imbalance in C/0 
atom numbers biases dust formation strongly toward either 
carbon-rich or oxygen-rich (including silicates). The ~ 1400 
K peak in AGN suggests an oxygen-rich environment from 
which the dust formed. 

The strength of our data set is that it allows us to ob- 



serve simultaneously the accretion disc and the hot dust 
emission. Therefore, we can derive for the first time mean- 
ingful covering factors for the dusty obscurer in AGN. If the 
ultraviolet radiation from the accretion disc emitted into 
the solid angle, defined by the dust distribution is com- 
pletely absorbed and re-emitted in the infrared, the dust 
covering factor is C = Q,/4m — J^^^L^dr^/ J^^^Li^du ~ 
0.4 ■ (/^Lhot/j^iacc), where i^Lhot and i^Lacc are the peak 
luminosities of the hot blackbody and accretion disc spec- 
trum, respectively (e.g.. lBarvainis|[r987l : iGranato fc Danesd 
1 19941 ). Note that whereas vLhot lies only slightly outside the 
observed spectral range and, therefore, is well-constrained 
by the data, the accretion disc peak luminosity is strongly 
model-dependent. In Fig. |9] we show the distribution of 
hot dust covering factors. We obtain values in the range 
of C ~ 0.01 - 0.6 and a mean of (C) = 0.07 ± 0.02. Our 
average value is a factor of ~ 6 lower than the average to- 



tal dust covering factor obtained by ISanders" et al.l (|1989H 
((C) = 0.40 ± 0.01), who considered the ratio between the 
entire integrated infrared luminosity (in the frequency range 
v ^ 10^^ — 10^"*'^ Hz) and the integrated accretion disc lumi- 
nosity for a large (~ 100 sources) sample of bright quasars. 

The relatively low hot dust covering factors could be a 
result introduced by the correction factors that we applied 
to the optical spectra (Section [3]), i.e., we have artificially 
raised the accretion disc fiux relative to the hot dust emis- 
sion. In order to address this concern, we have plotted in 
Fig. [To] the covering factors versus the optical correction 
factors. Whereas we see an envelope in the sense that the 
higher the optical correction factor, the lower the highest 
covering factor, overall we observe large hot dust covering 
factors (C > 0.1) in only a few cases (mostly in Ark 120 and 
NGC 4151) or, put differently, AGN states cluster around 
relatively low C values independent of the optical correction 
factor. 

4.3 Additional continuum components 

The observed AGN SEDs can be explained mainly by the 
sum of an accretion disc spectrum and emission from hot 
dust. However, two additional continuum components are 
expected to be present, which albeit weaker could alter 
the AGN continuum slope. The first is the diffuse contin- 
uum (thermal emission and scattering) from the BELR, 
which can be strong due to the relatively high gas densities 
found in these regions. The Balmer jump and continuum 
(A < 3646 A) were identified in AGN spectra long ago (e.g., 
iMalkan fc Sarge nt 1982) and predicted in e arly photoioniza- 
tion models as weU (e.g.. IWills et al.lll985h . If the BELR is 
emitting a Balmer continuum, it is also emitting a Paschen 
continuum (A < 8205 A), the strength of which could be 
important throughout the optical and near-IR spectrum. 
Other, weaker, recombination continua (Bracke t t, etc .) will 
also be present within the IR. iKorista fc Goad] ll200ll) used 
the photoionization code CLOUDY ( Ferland et al. Il998l ) and 
calculated the full diffuse BELR continuum assuming clouds 
with simple distributions in gas density and incident ioniz- 
ing flux. They found that after correcting for host galaxy 
light ~ 20% of the flux at 5100 A may be due to this com- 
ponent, which is dominated by the thermal continua (nearly 
40% at the Balmer jump; see Korista fc Goad for details). In 
a representation of the SED as in Fig. [5] (log i^/^ vs. log 
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the relative contribution of this component is expected to 
strengthen with increasing wavelength as the accretion disc 
spectrum falls, resulting in a flattening of the near-UV to 
near-IR AGN continuum spectrum as compared to an ac- 
cretion disc alone. 

A second additional component is due to the finite 
albedo of t he same dusty clouds em itting the thermal IR 
continuum. iKorista fc Ferlandl (| 19981) investigated the effec- 
tive albedo of a variety of gas clouds potentially present 
within the BELR. These contributions are included in the 
diffuse continuum just described. They also presented the ef- 
fective albedo of a high column density, dusty cloud with pa- 
rameters grossly approximating the clouds producing most 
of the hot dust emission. They found a ~ 20% albedo span- 
ning the wavelength range of A = 2500 A - 1 ^ra that de- 
clines significantly for wavelengths outside. While the de- 
tailed wavelength-dependent albedo depends on the grain 
size distributions and composition, their overall amplitudes 
should not differ greatly. Although the spatially-integrated 
scattered nuclear light will depend on factors such as the 
geometry of the scattering region and the observer's view- 
ing angle, contributions of ~ 10% to the AGN continuum 
are possible. Due to the finite grain albedo, the dusty torus 
will contribute additional light to the near-UV to near-IR 
spectrum that otherwise would not be expected from its 
thermal emission alone. As with the contribution from the 
diffuse continuum light from the BELR, this contribution 
will act to flatten (in log vfi, vs. log v) the AGN continuum. 

Due to the uncertainties in the absolute flux calibra- 
tion and the relative normalizations of our near-IR and op- 
tical spectra, we have not attempted to isolate the contri- 
butions of the diffuse BELR continuum and the dust scat- 
tering of nuclear light to the AGN continuum. However, do- 
ing so should provide additional physical and geometrical 
constraints on the BELR and inner dusty torus. Here we 
only point out that the strongest contributor to the diffuse 
BELR continuum, namely the Balmer jump and continuum 
are clearly visible in most of our optical spectra, forming 
together with Fc II the so-called 'small blue bump' (|Grandil 
[l982 : .Wills et al...l985. ) . 



5 THE AGN CONTINUUM VARIABILITY 

Both AGN continuum components sampled by our spec- 
tra are expected to be variable. The main cause of vari- 
ability will be a change in the accretion rate, which will 
lead to a change in the dust illumination. However, al- 
though the variability behaviour of the accretion disc has 
been studied extensively by optical reverberation cam- 
paigns, only few sources have bee n observed i n dust rever- 
boration programmes so far (e.g., 'Glass"l992'; 'Nelson"l996!; 
Oknyanskii ct al 199 9; Glass ,2004; Minozaki ct al. 200J; 
Suganuma et allbood : Ixoshida et alll20090 . Therefore, the 
exact location of the dusty obscurer as well as its geometry 
remain uncertain. In this section we discuss the first results 
of our coordinated optical and near-IR reverberation cam- 
paign. Based on these results, we will select the final sample 
for monitoring and develop a suitable observational strategy. 

We have multiple observation epochs available for 15/21 
sources listed in Table [5l of which five sources have more 
than two epochs. Defining the variability factor as the flux 
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Figure 11. Histogram of the variability factor measured at 1 
and calculated over a fixed period of one year. 
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Figure 12. Change in hot dust luminosity versus change in accre- 
tion disc luminosity over a fixed period of one year. The horizontal 
solid line indicates the locus of equal hot dust luminosity in both 
AGN states, with a deficit of hot dust observed in the high state 
for sources below this line. The diagonal solid line marks the locus 
of equality. Open and filled circles indicate sources with average 
hot dust radii greater and lower than one light year, respectively. 



change at 1 fim relative to the lowest-flux state in a fixed 
period of one year, we obtain the distribution shown in Fig. 
111! Note that this variability factor is different from that 
listed in Table [5] (column (6)), which refers to the total 
period between two epochs, and was calculated assuming 
uniform variability between the lowest-fiux epoch and the 
epoch with a higher fiux. The most variable of our sources 
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are 3C 273, HE 1228+013 and NGC 4151, which have at 
least one epoch with a one-year variability factor of > 2, fol- 
lowed by PG 0844-)-349, Ark 120, NGC 5548 and NG^ 4593, 
with variability factors of ~ 1.5 — 2. The source Mrk 817 
shows no variability and our three least variable sources are 
IRAS 1750-F508, Mrk 509 and H 2106-099. 

With our future near-IR reverberation programme we 
will be able to study also the variability behaviour of the hot 
dust and its dependence on the state of the ionising source. 
A first look gives Fig. 1121 where we have plotted the hot dust 
variability factor versus the accretion disc variability factor, 
with both values calculated for a fixed period of one year. As- 
suming that the hot dust absorbs and re-emits all the accre- 
tion disc luminosity, we have then estimated the average ra- 
dius of the hot dust component as -Rhot = \/ '^Lacc/^naT^^^, 
where a is the Stefan-Boltzmann constant, and have sepa- 
rated sources into those having -Rhot > 1 and < 1 light year 
(open and filled circles, respectively). 

Two important results become evident from Fig. 1121 
Firstly, for most sources we observe a response of the hot 
dust to the change in accretion disc luminosity, which is al- 
ways lower than the latter, as expected for a response lag. 
Most interestingly, however, is that a dust response is ob- 
served even in those sources with estimated hot dust radii 
well above one light year, i.e., for which we would expect 
no response. The only exception in this group is the source 
IRAS 1750-1-508. Secondly, for a few sources we observe a 
deficit of hot dust in the high state. This effect is most pro- 
nounced hi_the_sourceNG^4 for which it was observed 
also by iKoshida et all l|2009l ) using long-term optical and 
near-IR imaging and attributed to dust destruction. 

6 SUMMARY AND CONCLUSIONS 

We have used four epochs of quasi-simultaneous (within two 
months) near-IR and optical spectroscopy of 23 broad-line 
AGN to study the continuum spectral shape around 1 /im. 
Our main results can be summarized as follows. 

(i) The accretion disc spectrum appears to dominate the 
flux at ~ 1 /im, which allows us to derive a new relation that 
can be used to estimate AGN black hole masses. It is based 
on the near-IR virial product, defined here as the product 
between the width of the Pa/3 broad emission line and the in- 
tegrated 1 fim continuum luminosity. The dominance of the 
accretion disc spectrum at such long wavelengths means that 
the AGN state can be determined directly from the near-IR 
spectrum, making simultaneous optical spectroscopy for a 
reverberation programme unnecessary. 

(ii) An adequate subtraction of (in particular optical) 
host galaxy light reveals that the AGN continuum in the 
rest-frame frequency range of i/ ~ 10" - 10^^ Hz can be ap- 
proximated by the sum of mainly two emission components, 
a hot dust blackbody and an accretion disc spectrum. 

(iii) For the hot dust component we derive temperatures 
in the range of Thot '^-^ 1100 — 1700 K, which are typical 
values of the dust sublimation temperature, with a mean 
of (Thot) = 1365 ± 18 K. This mean value is close to the 
sublimation temperature of silicate dust grains, indicating 
that either carbonaceous dust is rare in AGN or that some 
other mechanism than dust sublimation sets the maximum 
dust temperature. The resulting hot dust covering factors 



are relatively low and in the range of C ~ 0.01 — 0.6, with a 
mean of (C) = 0.07 ± 0.02. 

(iv) Our preliminary variability studies have identified 
promising candidates for a future near-IR reverberation pro- 
gramme. Our three most variable sources in the near-IR are 
3C 273, HE 1228-^013 and NGC 4151. Furthermore, we have 
studied the response of the hot dust emission to changes in 
the accretion disc fiux. Most sources show the expected time 
lag, but a few sources have a deficit of hot dust in the high 
state, which indicates dust destruction. 

In our future work we will study the variability of the 
near-IR broad emission lines and constrain their physical 
conditions using detailed photoionisation models. In the 
longer term we plan to image our sample with current and 
future near-IR interferometers (|Elvis fc Karovskall2002l l. 
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APPENDIX A: ORIGINAL AND CORRECTED 
SPECTRAL ENERGY DISTRIBUTIONS 



This paper has been typeset from a T^^X/ M^jX file prepared 
by the author. 
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Figure Al. Rest-frame spectral energy distributions for the IRTF observing runs in 2004 May (black), 2006 January (green), 2006 June 
(red), and 2007 January (blue). Left and right panels show the original data and the data after host galaxy subtraction and spectral 
alignment were applied, respectively. See text for more details. 
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Figure Al - continued 



22 H. Landt et al. 




Figure Al - continued 
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Figure Al - continued 
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Figure Al - continued 
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Figure Al - continued 
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Figure Al - continued 
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Figure Al - continued 



